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We present a comprehensive black-hole event generator, BlackMax, which simulates the experi- 
mental signatures of microscopic and Planckian black-hole production and evolution at the LHC in 
the context of brane world models with low-scale quantum gravity. The generator is based on phe- 
nomenologically realistic models free of serious problems that plague low-scale gravity, thus offering 
more realistic predictions for hadron-hadron colliders. The generator includes all of the black-hole 
gray-body factors known to date and incorporates the effects of black-hole rotation, splitting be- 
tween the fermions, non-zero brane tension and black-hole recoil due to Hawking radiation (although 
not all simultaneously). 

The generator can be interfaced with Herwig and Pythia. 

The main code can be downloaded from 
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PACS numbers: 04.50.Gh, 04.70.Dy 
I. INTRODUCTION 

Models with TcV-scale quantum gravity 0, H, 0i Q 
offer very rich collider phenomenology. Most of them 
assume the existence of a three-plus-one-dimensional 
hypersurface, which is referred as "the brane," where 
Standard-Model particles are confined, while only gravity 
and possibly other particles that carry no gauge quantum 
numbers, such as right handed neutrinos can propagate 
in the full space, the so called "bulk" . Under certain as- 
sumptions, this setup allows the fundamental quantum- 
gravity energy scale, M,, to be close to the elcctroweak 
scale. The observed weakness of gravity compared to 
other forces on the brane (i.e. in the laboratory) is a 
consequence of the large volume of the bulk which di- 
lutes the strength of gravity. 

In the context of these models of TeV-scale quantum 
gravity, probably the most exciting new physics is the 
production of micro-black-holes in near-future accelera- 
tors like the Large Hadron Collider (LHC) Accord- 
ing to the "hoop conjecture" if the impact parame- 
ter of two colliding particles is less than two times the 
gravitational radius, r/j, corresponding to their center- 
of-mass energy (Ecu), a black- hole with a mass of the 
order of Eqm and horizon radius, r^, will form. Typi- 
cally, this gravitational radius is approximately Eqm/M^ 
Thus, when particles collide at center-of-mass energies 
above M^, , the probability of black- hole formation is high. 

Strictly speaking, there exist no complete calcula- 
tion (including radiation during the process of forma- 
tion and back-reaction) which proves that a black-hole 
really forms. It may happen that a true event horizon 



and singularity never forms, and that Hawking (or rather 
Hawking-like) radiation is never quite thermal. In Q this 
question was analyzed in detail from a point of view of an 
asymptotic observer, who is in the context of the LHC the 
most relevant observer. It was shown that though such 
observers never observe the formation of an event hori- 
zon even in the full quantum treatment, they do register 
pre-Hawking quantum radiation that takes away energy 
from a collapsing system. Pre-Hawking radiation is non- 
thermal and becomes thermal only in the limit when the 
horizon is formed. Since a collapsing system has only 
a finite amount of energy, it disappears before the hori- 
zon is seen to be formed. While these results have im- 
portant implications for theoretical issues like the infor- 
mation loss paradox, in a practical sense very little will 
change. The characteristic time for gravitational collapse 
in the context of collisions of particles at the LHC is very 
short. This implies that pre-Hawking radiation will be 
quickly experimentally indistinguishable from Hawking 
radiation calculated for a real black hole. Also, calcula- 
tions in Q indicate that the characteristic time in which 
a collapsing system losses all of its energy is very similar 
to a life time of a real black-hole. Thus, one may proceed 
with a standard theory of black-holes. 

Once a black-hole is formed, it is believed to decay via 
Hawking radiation. This Hawking radiation will consist 
of two parts: radiation of Standard-Model particles into 
the brane and radiation of gravitons and any other bulk 
modes into the bulk. The relative probability for the 
emission of each particle type is given by the gray-body 
factor for that mode. This gray-body factor depends on 
the properties of the particle (charge, spin, mass, momen- 
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turn), of the black- hole (mass, spin, charge) and, in the 
context of TeV-scale quantum gravity, on environmental 
properties - the number of extra dimensions, the location 
of the black- hole relative to the brane (or branes), etc.. In 
order to properly describe the experimental signatures of 
black-hole production and decay one must therefore cal- 
culate the gray-body factors for all of the relevant degrees 
of freedom. 

There are several black-hole event generators available 
in the literature Q based on particular, simplified mod- 
els of low-scale quantum gravity and incorporating lim- 
ited aspects of micro-black-hole physics. Unfortunately, 
low-scale gravity is plagued with many phenomenolog- 
ical challenges like fast proton decay, large nfi oscilla- 
tions, flavor-changing neutral currents and large mixing 



between leptons |ld . Illj . For a realistic understanding of 



the experimental signature of black hole production and 
decay, one needs calculations based on phenomenologi- 
cally viable gravity models, and incorporating all neces- 
sary aspects of the production and evolution of the black- 
holes. 

One low-scale gravity model in which the above men- 
tioned phenomenological challenges can be addressed is 



the split-fermion model In this model, the Stan- 

dard Model fields are confined to a "thick brane" , much 
thicker than M~^. Within this thick brane, quarks and 
leptons are stuck on different three-dimensional slices 
(or on different branes), which are separated by much 
more than M~^. This separation causes an exponen- 
tial suppression of all direct couplings between quarks 
and leptons, because of exponentially small overlaps be- 
tween their wave-functions. The proton decay rate will 
be safely suppressed if the spatial separation between 
quarks and leptons is greater by a factor of at least 10 
than the widths of their wave functions. Since AB = 2 
processes, like nn oscillations, are mediated by operators 
of the type uddudd, suppressing them requires a further 
splitting between up-type and down-type quarks. Since 
the experimental limits on AB = 2 operators are much 
less stringent than those on AB — 1 operators, the u and 
d-type quarks need only be separated by a few times the 
width of their wave functions [l2| ■ 

Current black-hole generators assume that the black- 
holes that are formed are Schwarzschild-like. However, 
most of the black-holes that would be formed at the LHC 
would be highly rotating, due to the non-zero impact pa- 
rameter of the colliding partons. Due to the existence 
of an ergosphere (a region between the infinite redshift 
surface and the event horizon) , a rotating black-hole ex- 
hibits super-radiance: some modes of radiation get am- 
plified compared to others. The effect of super-radiance 
[l3| is strongly spin-dependent, with emission of higher- 
spin particles strongly favored. In particular the emission 
of gravitons is enhanced over lower-spin Standard-Model 
particles. Since graviton emission appears in detectors as 
missing energy, the effects of black-hole rotation cannot 



be ignored. Similarly, black-holes may be formed with 
non-zero gauge charge, or acquire charge during their 
decay. This again may alter the decay properties of the 
black-hole and should be included. 

Another effect neglected in other generators is the re- 
coil of the black-hole. A small black-hole attached to a 
brane in a higher-dimensional space emitting quanta into 
the bulk could leave the brane as a result of a recoil^ . In 
this case, visible black-hole radiation would cease. Alter- 
nately, in a split-brane model, as a black-hole traverses 
the thick brane the Standard-Model particles that it is 
able to emit will change depending on which fermionic 
branes are nearby. 

It is also the case that virtually all the work in this field 
has been done for the idealized case where the brane ten- 
sion is negligible. However, one generically expects the 
brane tension to be of the order of the fundamental en- 
ergy scale, being determined by the vacuum energy con- 
tributions of brane-localized matter fields [ij]. As shown 



15[, finite brane tension modifies the standard gray- 



body factors. 

Finally, it has been suggested [3l that more common 
than the formation and evaporation of black-holes will be 
gravitational scattering of parton pairs into a two-body 
final state. We include this possibility. 

Here we present a comprehensive black- hole event gen- 
erator, BlackMax, that takes into account practically all 
of the above mentioned issues^, and includes almost all 
the necessary gray-body factors'^. Preliminary studies 
show how the signatures of black-hole production and 
decay change when one includes splitting between the 
fermions, black-hole rotation, positive brane tension and 
black-hole recoil. Future papers will explore the implica- 
tions of these changes in greater detail. 

In sectioniniand lllll we discuss the production of black- 
holes and the gray-body factors respectively. The evap- 
oration process and final burst of the black-holes is dis- 
cussed in section Hvl and FVl Sections I VII and describe the 
input and output of the generator. Section IVIII shows 
some characteristic distributions of black-holes for differ- 
ent extra dimension scenarios. The reference list is exten- 
sive, refiectin g th e great interest in the 
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but by no means complete. 



^ Although if the black-hole earries gauge charge it will be pre- 
vented from leaving the brane. 
^ although not necessarily simultaneously 

^ Except in the one case of the graviton gray-body factor for a 
rotating black-hole, where the calculation has yet to be achieved. 
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II. BLACK-HOLE PRODUCTION 

We assume that the fundamental quantum-gravity en- 
ergy scale A'/* is not too far above the electroweak scale. 
Consider two particles colliding with a center-of-mass en- 
ergy EcM ■ They will also have an angular momentum J 
in their center-of-mass (CM) frame. By the hoop conjec- 
ture, if the impact parameter, 6, between the two collid- 
ing particles is smaller than the diameter of the horizon 
of a ((i + l)-dimensional black-hole (where d is the total 
number of space-like dimensions) of mass M — Ecm and 
angular momentum J, 



b < 2rhid,M, J), 



(1) 



then a black-hole with will form. The cross section 
for this process is approximately equal to the interaction 
area 7r(2r;i)^. 

In Boyer-Lindquist coordinates, the metric for a (d-l-l)- 
dimensional rotating black-hole (with angular momen- 
tum parallel to the uj in the rest frame of the black-hole) 
is: 



ds^ = 1 - 



4-d 



S(r,( 



■ 2 n 2 , 2 

sm [ r + a 



4-d 



S(r,( 



2a sin2 e^^dtdS - ^^dr^ 



S](r, 6l)d6l^ - r^cos^Od^-^n 



(2) 



where /i is a parameter related to mass of the black-hole, 
while 



and 



A = r^ + a^ -fir^-". 



The mass of the black-hole is 

id-l)Ad-i 



M 



and 



J 



2Ma 



d-1 

is its angular momentum. Here, 

Ad-i 



r(d/2) 



(3) 
(4) 

(5) 
(6) 

(7) 



is the hyper-surface area of a (d — l)-dimensional unit 
sphere. The higher-dimensional gravitational constant 
Gd is defined as 



Gd — 



4Mt 



(8) 
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FIG. 1: Horizon radius (in GeV"^) of a non-rotating black- 
hole as a function of mass for 4-10 spatial dimensions. 



The horizon occurs when A = 0. That is at a radius 
given implicitly by 



Here 



^(d) ^ ^l/id-2) 



(9) 



(10) 



is the Schwarzschild radius of a (d-l-l)-dimensional black- 
hole, i.e. the horizon radius of a non- rotating black- hole. 
Equation 1101 can be rewritten as: 



ll/(d-2) 



where 



k{d) 



2rf-3^(d-6)/2 



r[d/2] 



-|l/(d-2) 



1 



(11) 



(12) 



Figure [T] shows the horizon radius as a function of 
black-hole mass for d from 4 to 10. We see that the hori- 
zon radius increases with mass; it also increases with d. 
Figure [H shows the Hawking temperature of a black-hole 



d-2 
47rr,i 



(13) 



as a function of the black-hole mass for d from 4 to 10. 
The Hawking temperature is a measure of the character- 
istic energies of the particles emitted by the black-hole. 
Th decreases with increasing mass. However, the be- 
haviour of Th with changing d is complicated, reflecting 
the competing effect of an increasing horizon radius and 
an increasing d — 2 in equation 1131 

For the model with with non-zero tension brane, the 
radius of the black-hole is defined as 



„(*) 



(14) 
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FIG. 2: Hawking temperature(in GeV) of a non-rotating 
black-hole as a function of mass for 4-10 spatial dimensions. 



FIG. 4: Hawking temperature(in GeV) of a black-hole as a 
function of mass for different B in d=5 spatial dimensions. 
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FIG. 3: Horizon radius (in GeV ^) of a black- hole as a func- 
tion of mass for different B in d=5 spatial dimensions. 



FIG. 5: Horizon radius (in GeV~^) of a rotating black-hole 
as a function of mass for different angular momentum in d=5 
spatial dimensions. Angular momentum J is in unit of h. 



with B the deficit-angle parameter which is inverse pro- 
portional to the tension of the brane. 

Figure [3] shows the horizon radius function of 

black-hole mass for the model with non-zero tension 
brane. As the deficit-angle parameter increases, the size 
of the black-hole increases. 

Figure [4] shows the Hawking temperature of a black- 
hole for the model with non-zero tension brane. The 
Hawking temperature decreases as the deficit angle de- 
creases. 

Figure [5] shows the horizon radius as a function of 
black-hole mass for a rotating black-hole. The angular 
momentum decreases the size of the horizon and increases 
the Hawking temperature (see figures [5] and [6|) . 

If two highly relativistic particles collide with center-of- 
mass energy Ecm, and impact parameter b, then their 
angular momentum in the center-of-mass frame before 
the collision is Lin = bEcAi/'^- Suppose for now that 
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FIG. 6: Hawking temperature(in GeV) of a rotating black- 
hole as a function of mass for different angular momentum in 
d=5 spatial dimensions. Angular momentum J is in unit of 

h. 
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the black-hole that is formed retains all this energy and 
angular momentum. Then the mass and angular momen- 
tum of the black-hole will be Mi„ = Ecu and Jin = L.^. 
A black-hole will form if: 



b < Ecm/2) . 



(15) 



We see that bmax is a function of both Ecm and the 
number of extra dimensions. 

We can rewrite condition (fTSl) as 



[Ecm ; d) 



CM) 



1 



V 2 ) 



(16) 



There is one exception to this condition. In the case 
where we are including the effects of the brane tension, 
the metric (and hence gray-body factors) for a rotating 
black-hole are not known. In this case we consider only 
non-rotating black-holes. Therefore, for branes with ten- 
sion 



"max [Ecm, a) 



2ri''\Ecu). 



(17) 



Also, for branes with tension only the d = 5 metric is 
known. 

At the LHC, each proton will have E = 7TeV in the 
CM frame. Therefore, the total proton-proton center of 
mass energy will be ^/s = 14TcV. However, it is not the 
protons that collide to make the black-holes, but the par- 
tons of which the protons are made. If two partons have 
energy vE and much greater than their respective 
masses, then the parton-parton collision will have 

s' = bi+Pjf = \v{E,E) + -{E,-E)\^ ^4uE^ ^us. 

(18) 

We define a quantity Q' 



Q' = Ecm = Vs' = Vus 



(19) 



The centcr-of-mass energy for the two colliding partons 
will be y/us, as will be the 4- momentum transfer Q'^. 
The largest impact parameter between the two partons 
that can form a black-hole with this mass will therefore 
be bmax{,\fus] d) , as given by equation [TBI 

The total proton-proton cross section for black-hole 
production is therefore 



pp^BH 



;d, A4)= / du — TT \bmaxiVus; d)]' 

J Mi Is Ju V ^ 



(20) 



Here fi{v,Q') is the i-th parton distribution function. 
Loosely this is the expected number of partons of type i 
and momentum vE to be found in the proton in a colli- 
sion at momentum transfer Q' . 



In [l^ it is argued that strong gravity effects at ener- 
gies close to the Planck scale will lead to an increase in 
the 2^2 cross section via the exchange of Planckian 
"black-holes" (by which any quantum gravity effect or 
resonance is meant). Final states with high multiplicities 
are predicted to be suppressed. Although the intermedi- 
ate state is created in the strong gravity regime, it is not 
a conventional microscopic black hole. The state is not 
stable. Thermal Hawking radiation does not take place. 
Especially since inelastic collisions increase the energy 
loss, the threshold for creating stable black-holes shifts 
to even higher values. Thus 2^2 scattering may be 
the most important signal in the LHC instead of black- 
holes evaporating via Hawking radiation. One should 
find that the cross section for two-body final states sud- 
denly jumps to a larger value, as the energy reaches the 
quantum-gravity scale. We calculate the cross section of 
two-body final states by replacing Trfe^^j^. in equation ([20]) 
with 



T^h„iax{y S' > Mmin) 



where 



P2 

<N> 
P 



'<N> 



< N >' 



E 

ATrkjd) MBH .id-i)/{d-2) 

Ec^g.r.C(3)r(3) 
Ecj,a>.c(4)r(4) 

1 f ai{ijj)iLp'dw 



47rr2 ./ e'^/^ ± 1 
1 f ai{ijj)uj^duj 



oj du) 
uj^duj , 



(21) 

(22) 

(23) 
(24) 
(25) 
(26) 



Here Ci is the number of internal degrees of freedom 
of particle species i, = 1 and fi = 1 for bosons, and 
= 3/4 and /, = 7/8 for fcrmions 

Figure[7]shows the cross section for non-rotating black- 
holes on a tensionless brane as a function of mass for dif- 
ferent numbers of spatial dimensions. The cross section 
increases with the number of spatial dimensions. 

Figure [8] shows the cross section for non-rotating black 
holes on a brane with positive tension as a function 
of mass for various deficit angle parameter, B. The 
cross section increases as the tension increases (as B de- 
creases). 

Figure [9] shows the cross section for non- rotating black 
holes as function of the number of split-fermion space 
dimensions, Us- When a pair of partons are separated in 
the extra-dimensions they must approach more closely 
in the ordinary dimensions in order to form a black-hole. 
Thus the effective cross-section for black-hole formation 
in collisions is decreased. This effect become more severe 
as Tis increases because the partons are more likely to be 



6 



^0■ : 



10' 



" 10 : 



10 ' 



M =1TeV 

M =3TeV 

. .M=5TeV 
M =5TeV M =14TeV 
Single brane model 



5 e 7 8 9 10 
Number of spatial dimensions 



10' 



10' -- 



M =1TeV 

M =3TeV 

■ M =5TeV 
M,,„=5TeV M, =MTey 
d=10 

Splitfernnionmodel 



2 3 4 5 
Number of splitting space 



FIG. 7: Cross section for production of a black-hole (rotating 
or non-rotating) as a function of the number of spatial dimen- 
sions, for a tensionless brane, with no fermion brane-splitting. 



FIG. 9: Cross section for production of a non-rotating black- 
hole as a function of the number of fermion brane-splitting 
dimensions for d = 10. 
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FIG. 8: Cross section for production of a non-rotating black- 
hole as a function of the deficit angle parameter for d = 5 and 
Us = 2. 



more widely separated in the extra dimensions, therefore 
the cross section decreases with increasing Us . 

Figure (TU] shows the cross section as a function of the 
chosen minimum black-hole mass. The parton distri- 
bution functions strongly suppress the events with high 
black-hole masses. 

Figure [Tl] shows the cross section for the two-body 
final-state scenario as a function of the number of spatial 
dimensions, for Mmin = = ITeV, M„ii„ = = 
3TeV and Mmin = M^, = 5TeV. It increases with the 
number of spatial dimensions. 



A. Black-Hole Formation in BlackMax 



FIG. 10: Cross section for formation of a black-hole (rotating 
or non-rotating) as function of the minimum mass of black- 
hole, for a zero-tension brane, with no fermion brane-splitting. 
The vertical lines are the error bars. 

protons of center-of-mass-energy ^/s, is given by 

P{Q')- f'-J2Mv,Q')M-,Q')- (27) 

According to the theory, there will be some minimum 
mass for a black-hole. We expect Mmin ^ M^,, but leave 
Mmin > M^, as a free parameter. Therefore, a black-hole 
will only form if m > (Mmin/Q)^- The type of partons 
from which a black-hole is formed determines the gauge 
charges of the black-hole. Clearly, the probability to cre- 
ate a black-hole in the collision of any two particular 
partons i and j with energies (momenta) vE and is 
given by 



Within BlackMax, the probability of creating a black- 
hole of center-of-mass energy y^ws, in the collision of two 



PivE,—,i,j) = Mv,Q')fj{u/v,Q') (28) 



7 





10 ■ 






10'- 




M =r,l = IT*-.- 

- - M =M =3 7*^: 








. .r,i=M =frT«v 


D 






n.=0 H ,|^_=14T*V 




10' - 
















10'- 


_. ■-■ 






10 ■ 








10 '- 


1 I 1 ■ 1 ■ 





6 7 S 9 10 

Number of spatial dimensions 



FIG. 11: Cross section for the two-body final-state scenario 
as a function of number of spatial dimensions where Mmin = 
M, = ITeV, M„ii„ = M. = 3TeV and M„i„ = 5TeV. 



The energies and types of the two colhding partons de- 
termine their momenta and affect their locations within 
the ordinary and extra dimensions. For protons moving 
in the z-direction, we arbitrarily put one of the partons at 
the origin and locate the second parton randomly within 
a disk in the xy-plane of radius bmax{EcM = \/us] d). 

We must, however, also take into account that the par- 
tons will be separated in the extra dimensions as well. 
Each parton type is given a wave function in the ex- 
tra dimension. For fermions, these wave functions are 
parametrized by their centers and widths which are in- 
put parameters (cf. Fig [T4)) . In the split-fermion case, the 
centers of these wave functions may be widely separated; 
but even in the non-split case, the wave functions have 
non-zero widths. For gauge bosons, the wave functions 
are taken to be constant across the (thick) brane. 

The output from the generator (described in greater 
detail below) includes the energies, momenta, and types 
of partons that yielded black- holes. The locations in time 
and space of the black- holes are also output. 

The formation of the black-hole is a very non-linear 
and complicated process. We assume that, before set- 
tling down to a stationary phase, a black-hole loses some 
fraction of its energy, linear and angular momentum. We 
parameterize these losses by three parameters: 1 — Je, 
1 ~ fp and 1 — /l- Thus the black- hole initial state that 
we actually evolve is characterized by 



and angular momentum that are retained by the station- 
ary black-hole. We expect that most of the energy lost in 
the non-linear regime is radiated in the form of gravita- 
tional waves and thus represents missing energy. Yoshino 
and Rychkov Q have calculated the energy losses by 
numerical simulation of collisions. Their results will be 
incorporated in a future upgrade of BlackMax. 

For a small black-hole, the numerical value for the an- 
gular momentum is of the order of several fi. In that 
range of values, angular momentum is quantized. There- 
fore a black-hole cannot have arbitrary values of angu- 
lar momentum. We keep the actual angular momentum 
of the black-hole, J, to be the nearest half-integer, i.e. 
2J=[2J'-t-i]. 

The loss of the initial angular momentum in the non- 
linear regime has as a consequence that the black-hole 
angular momentum is no longer in the transverse plane 
of the colliding protons. We therefore introduce a tilt in 
the angular-momentum 

e = cos-\^ ). (30) 

Figure [T2l illustrates this geometry. 

In this version of the generator, we have assumed that 
the angular-momentum quantum numbers of the black- 
hole were ( J, Jm = J)'^- We next randomly choose an 
angle 0, and then reset the angular-momentum axis to 




E 

P. 
J' 



Ein Ie ; 

Ez in fp 5 



(29) 



FIG. 12: Angular momentum tilt geometry. 



where Ein , Pz in and Lin are initial energy, momentum 
and angular momentum of colliding partons, while Je, fp 
and /l are the fractions of the initial energy, momentum 



Future versions of the generator may randomize the choice of 

J m • 
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III. GRAY-BODY FACTORS 



Once the black-hole settles down to its stationary con- 
figuration, it is expected to emit semi-classical Hawk- 
ing radiation. The emission spectra of different parti- 
cles from a given black-hole depend in principle on the 
mass, spin and charge of the black-hole, on the "en- 
vironment"^ and on the mass and spin of the partic- 
ular particle. Wherever possible we have made use of 
the correct emission spectrum often phrased in terms of 
the gray-body factor for black- holes in 3+1-dimensional 
space-time. In most cases, these were extant in the lit- 
erature, but we have calculated the spectra for the split- 
fermion model ourselves, and reproduced existing spectra 
independently. The sources of the gray-body factors are 
summarized in Table HI 

• Non-rotating black-hole on a tensionless hrane: 
For a non-rotating black-hole, we used previously 
known gray-body factors for spin 0, 1/2 and 1 fields 
in the brane, and for spin 2 fields (i.e. gravitons) 
in the bulk. 

• Rotating black-hole on a tensionless brane: For ro- 
tating black- holes, we used known gray-body fac- 
tors for spin 0, 1/2 and 1 fields on the brane. The 
correct emission spectrum for spin 2 bulk fields 
is not yet known for rotating black-holes, we cur- 
rently do not allow for the emission of bulk gravi- 
tons from rotating black-holes. As discussed be- 
low, this remains a serious shortcoming of cur- 
rent micro-black-hole phenomenology, since super- 
radiance might be expected to significantly increase 
graviton emission from rotating black-holes, and 
thus increase the missing energy in a detector. 

• Non-rotating black-holes on a tensionless brane 
with fermion brane splitting: In the split-fermion 
models, gauge fields can propagate through the 
bulk as well as on the brane, so we have calculated 
gray-body factors for spin and 1 fields propagat- 
ing through the bulk, but only for a non-rotating 
black-hole for the split-fermion model. These are 
shown in Figures [50ll611 

• Non-rotating black-holes on a non-zero tension 
brane: The bulk gray-body factors for a brane 
with non-zero tension are affected by non-zero ten- 
sion because of the modified bulk geometry (deficit 
angle). We have calculated gray-body factors for 
spin 0, 1 and 2 fields propagating through the bulk. 



dimensionality and geometry of the bulk, brane tension, location 
of fermionic branes 



again only for the non-rotating black-hole for a 
brane with non-zero tension and d = 5. 

• Two particle final states: We use the same gray- 
body factors as a non-rotating black-hole to cal- 
culate the cross section of two-particle final states 
(excluding gravitons) . 

In all cases, the relevant emission spectra are loaded 
into a data base as described in appendices A. 

IV. BLACK-HOLE EVOLUTION 



The Hawking radiation spectra are calculated for the 
black-hole at rest in the center-of-mass frame of the col- 
liding partons. The spectra are then transformed to the 
laboratory frame as needed. In all cases we have not 
(yet) taken the charge of the black-hole into account 
in calculating the emission spectrum, but have included 
phenomenological factors to account for it as explained 
below. 

The degrees of freedom of the Standard-Model parti- 
cles are given in Table |TT1 Using the calculated Hawking 
spectrum and the number of degrees of freedom per par- 
ticle, we determine the expected radiated fiux of each 
type of particle as a function of black-hole and environ- 
mental properties. For each particle type i we assign to 
it a specific energy, hcoi with a probability determined by 
that particle's emission spectrum. (The particle "types" 
are listed in Table [Til) 

Assume a black-hole with mass M^fi emits a massless 
particle with energy huji. The remaining black-hole will 
have energy and momentum like 

{Mbh - f^i,-fi^i) (31) 

Here we ignore the other dimensions. We use a classical 
model to simulate the events. The mass of the remaining 
black-hole should remain positive. So from equation (|3ip 

one gets 

huj, < Mbh/2. (32) 

Combining this with the observation that energy of a 
particle is larger than its mass, leads us to require that 

M^ < hwi (33) 

We next need to determine whether that particle with 
that energy is actually emitted within one generator 
time-step At. The time-step itself is an input param- 
eter (cf. figure [T4|) . We choose a random number Nr 
from the interval [0, 1]. Given Lpi, the total number fiux 
of particles of type i, and A'^;, the number of degrees of 
freedom of that particle type, the particle will be emitted 
if 

LptN.At > Nr. (34) 
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FIG. 13: The emitted fermion intensity (normalized to one) 
as a function of the distance between the black-hole and the 
center of the gaussian distribution of a fermionic brane. As a 
black-hole increases the distance from a fermionic brane due 
to recoil the intensity of the emitted fermions of that type falls 
down quickly. The radius of black-hole is set to be 2M^^ . The 
width of the fermionic brane is . The plot is for the case 
of one extra dimension. 

In the single-brane model, Lpi is derived from the 
power spectrum of the Hawking radiation. In the 
split-fermion model, we include a suppression factor for 
fermions. The factor depends on the overlap between 
the particular fermion brane and the black-hole when 
the black-hole is not located on that fermion's brane. 
Fig. [13] shows how the spectrum of emitted particles 
changes as the black-hole drifts away from the center of 
the fermion brane. As a black-hole increases its distance 
from a fermion brane due to recoil, the intensity of the 
emitted fermions of that type declines quickly. 

If the particle is to be emitted, we choose its angular- 
momentum quantum numbers (/, m) according to: 

P,^{i,l,m,E)^Lu^^{E)/ Lu>.r,AE)- (35) 

I' ,m' 

Here Pem(i, ^, m, i?) is the probability that a type i 
particle with quantum numbers {l,m) will be emitted. 
Lij^m{E) is the emission spectrum of a particle of type 
i with quantum numbers (l,m). This step is omitted in 
the case of non-rotating black-holes since we do not fol- 
low the angular-momentum evolution of the black-hole. 

Once the quantum numbers of the emitted particle arc 
determined, we calculate the direction of emission ac- 
cording to the corresponding spheroidal wave function: 

Pem(e,0) - \^imiO,(b)\^smeAeA(b (36) 

Here Pe,„(^, is the probability of emission in the 
{0, (j)) direction for the angular quantum numbers (£, m). 
^im{(^,4') is the (properly normalized) spheroidal wave 



function of the mode with those angular quantum num- 
bers. 

Once the energy and angular-momentum quantum 
numbers are determined for the i-th particle type, then, if 
that particle type carries SU(3) color we assign the color 
randomly. The color is treated as a three-dimensional 
vector Ci = (r, 5, 17) , in which a quark's color- vector is ei- 
ther (1, 0, 0), (0, 1, 0), or (0, 0, 1). Similarly, an anti-quark 
has —1 entries in its color- vector. A gluon's color- vector 
has one -1-1 entry, and one —1 entry. 

The black-hole gray-body factors which quantify the 
relative emission probabilities of particles with different 
spin are calculated for a fixed background, i.e. assuming 
that the black-hole metric does not change during the 
emission process. However, as the black-hole emits par- 
ticles, the spin and charge of the black-hole do change. 

A. Electric and Color Charge Suppression 

A charged and highly rotating black-hole will tend to 
shed its charge and angular momentum. Thus, emis- 
sion of particles with charges of the same sign as that 
of the black-hole and angular momentum parallel to the 
black-hole's will be preferred. Emission of particles that 
increase the black-hole's charge or angular momentum 
should be suppressed. The precise calculation of these 
effects has not as yet been accomplished. Therefore, to 
account for these effects we allow optional phenomeno- 
logical suppression factors for both charge and angular 
momentum. 

The following charge-suppression factors can currently 
be used by setting parameter 19 (cf. section IVl|) equal to 
2. 

F« = exp(CQg'''Q'=") (37) 
Fl = exp(C3cfcr) a = r,b,g. (38) 

Q^^ is the electromagnetic charge of the black-hole, Q*^™ 
is the charge of the emitted particle; c^'', is the color 
value for the color a, with a = r, b, g, of the black-hole, 
and c^™, is the color value for the color a, with a = r, b, g, 
of the emitted particle. Cg and C3 are phenomenological 
suppression parameters that are set as input parameters 
of the generator. 

We estimate = 0{aem) and ^3 = 0(as), where aem 
and as are the values of the electromagnetic and strong 
couplings at the Hawking temperature of the black-hole. 
Note that we currently neglect the possible restoration 
of the electroweak symmetry in the vicinity of the black- 
hole when its Hawking temperature is above the elec- 
troweak scale. Clearly, since Uem — 10~^ we do not 
expect electromagnetic (or more correctly) electroweak 
charge suppression to be a significant effect. However, 
since as(l TeV) ~ 0.1, color suppression may well play a 
role in the evolution of the black-hole. 
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Once we have determined the type of particle to be 
emitted by the black-hole, we draw a random number 
Nj. between and 1 from a uniform distribution. If 
Nj. > F'^ then the emission process is allowed to occur, 
if Nr < F'^ then the emission process is aborted. We re- 
peat the same procedure for color suppression factor, F^. 
Thus, particle emission which decreases the magnitude 
of the charge or color of the black-hole is unsuppressed; 
this suppression prevents the black-hole from acquiring a 
large charge/color, and gives preference to particle emis- 
sion which reduces the charge/color of the black- hole. 



B. Movement of the Black-Hole during 
Evaporation 



of the black-hole. However, as we shall see, the black- 
hole never becomes highly relativistic, so the recoil does 
not significantly affect the shape of the spectrum. 

As the lifetime of a small black-hole is relatively short, 
and its recoil velocity non-relativistic, it does not move 
far from its point of creation. However, even a recoil of 
the order of one fundamental length ^ in the bulk 

direction could dislocate the black-hole from the brane^. 
In single-brane models this would result in apparent miss- 
ing energy for an observer located on the brane able 
to detect only Standard Model particles. In the split- 
fermion model, as the black-hole moves off or on partic- 
ular fcrmion branes, the decay channels open to it will 
change. 



We choose the direction of the momentum of the emit- 
ted particle (Pg) according to equation ([36]) in the center- 
of-mass frame and then transform the energy and mo- 
mentum to their laboratory frame values hw' and P'^. The 
black-hole properties (energy, momentum, mass, colors, 
and charge) arc then accordingly updated for the next 
time step: 



E{t + At) 
P{t + M) 



m 
Pit) 



huj' 
P' 



M{t + At) 
*^{t + At) 
t^{t + At) 



= ^ E{t + At)^ - P{t- 



^bh 



Q'''{t)-Q^ 



(39) 

(40) 

At)^ (41) 
(42) 
(43) 



Here c*'' is the color 3-vector of the black-hole and can 
have arbitrary integer entries. 

Due to the recoil from the emitted particle, the black- 
hole will acquire a velocity v and move to a position x: 



v{t) 
x{t + At) 



x{t) +v{t)At 



(44) 
(45) 



Since fermions are constrained to live on the 3-1-1- 
dimensional regular brane, the recoil from fermions is 
not important. Only the emission of vector fields, scalar 
fields and gravitons gives a black hole momentum in ex- 
tradimension. Once a black hole gains momentum in 
extradimension, it is able to leave the regular brane if 
it carries no gauge charge. In the split fermion case, it 
can move within the mini-bulk even if it carries gauge 
charge. In the case of rotating black holes, because the 
gray-body factor for gravitons is not yet known, gravi- 
ton emission is turned off in the generator and the black 
holes experience no bulk recoil. 

Recoil can in principle change the radiation spectrum 
of the black- hole in two ways. First, the spectrum will 
not be perfectly thermal or spherically symmetric in the 
laboratory frame, but rather boosted due to the motion 



C. Rotation 

Since two colliding particles always define a single 
plane of rotation, rotating black-holes are formed with 
a single rotational parameter. For two particles colliding 
along the z-axis, there should be only one rotation axis 
perpendicular to the z-axis. However, due to angular- 
momentum loss both in the formation process, and sub- 
sequently in the black-hole decay, three things can hap- 
pen: i) the amount of rotation can change, ii) the rota- 
tion axis can be altered, and iii) more rotation axes can 
emerge, because there are more than three spatial dimen- 
sions. Also, if the colliding particles have a non-zero im- 
pact parameter in bulk directions^ the plane of rotation 
will not lie entirely in the brane direction. Because solu- 
tions do not exist for rotating black-holes with more than 
one rotation axis, we forbid the emergence of secondary 
rotation axes. We do, on the other hand, allow the sin- 
gle rotation axis of the black-hole to evolve. However, no 
gray-body factors are known if the single rotation axis ac- 
quires components in the extra dimensions, therefore we 
limit the rotation axis to the brane dimensions. Relaxing 
these limitations is a subject for future research. 

We next must determine the rotational axis of the 
black-hole. The rotation parameter of a black-hole with 
angular- momentum quantum numbers {j,j) is taken to 
be 



J n + 2 
M 2 ' 



(46) 



where J = \/ j{j + 1)^- The direction of the black-hole 
angular momentum is taken to be 



(47) 



^ This is very unlikely because most of the black-holes have gauge 
charges. 

Due to the finite thickness of the single brane or splitting between 
the quark branes. 
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where l± is a unit vector in the plane perpendicular to 
u!. We chose the direction of /j^ randomly. 

When the black-hole emits a particle with angular- 
momentum quantum numbers (/,m), there are several 
possible final states in which the black-hole can end up. 
We use Clebsch-Gordan coefficients to find the probabil- 
ity of each state. 



|j, j>=Xi^O',j;^m,j',i- m)\l,m> m > (48) 

j'<3+l 

We use |C(j, j;l,m, j', j — m)\'^ as the probability that 
the new angular-momentum quantum numbers of the 
black-hole will be — m). From angular-momentum 
conservation J = L + J' , we can calculate the tilt angle 
of the black-hole rotation axis as: 

cos0^i^l±^}±M±2l=I^. (49) 

We randomly choose a direction with the tilt angle 9 
as a new rotation axis and change quantum numbers to 
(/,/)• 

In calculating the gray-body factors, the black hole is 
always treated as a fixed unchanging background. The 
power spectrum of emitted particles can be calculated 
from 



dE 
'dt 



dto 



exp{{ui — mVL)/TH) =F 1 27r 



(50) 



Here I and m are angular momentum quantum numbers. 
Lo is the energy of the emitted particle. D, is definied by 



n = 



{l + al)rh 



(51) 



The exponential factor in the denominator of (jSP]) 
causes the black hole to prefer to emit high angular mo- 
mentum particles. However, since the TeV black holes 
are quantum black holes, the gray-body factors should 
really depend on both the initial and final black-hole pa- 
rameters. The calculation of the gray-body spectra on 
a fixed background can cause some problems. In par- 
ticular, in the current case, the angular momentum of 
the emitted particle (as indeed the energy) may well be 
comparable to that of the black hole itself. There should 
be a suppression of particle emission processes in which 
the black hole final state is very different from the initial 
state. We therefore introduce a new phenomenological 
suppression factor, parameter 17, to reduce the probabil- 
ity of emission events in which the angular momentum 
of the black hole changes by a large amount. 

If parameter 17 is equal to 1 (cf. section |Vl|, we do 
not take into account the suppression of decays which 
increase the angular momentum of the black-hole. If we 



are using A Area suppression (parameter 17 equal to 2) 
then 



exp(a(rr(t + At)7rr(t)^-l)) 



(52) 



If we are using J^/i suppression (parameter 17 equal to 3) 
then 



F^=exp(-a|J"'(i + AOI) 



(53) 



If we are using AJb^, suppression (parameter 17 equal to 
4) then 

F"^ = expi^ai-J'^t + At) - J^^m. (54) 

We might expect Cl ^ 1, however there is no detailed 
theory to support this; as indeed there is no detailed 
theory to choose among these three phenomenological 
suppression factors. It is also worth noting that, while for 
d ~ 3 and d ~ A there is a maximum angular momentum 
that a black-hole of a given mass can carry, for d > 5 
there is no such upper limit. We do impose that a < Rs/2 
for d = 3 and a < Rg for d ~ 4. 

As for the charge and color suppression, we choose a 
random number Nr between and 1. If Nr > then 
the particle emission is aborted. 

The procedure described in this section is then re- 
peated at each time step with each particle type, and 
then successive time steps are taken until the mass of 
the black- hole falls below M*. In practice, the time step 
should be set short enough that in a given time step the 
probability that particles of more than one type are emit- 
ted is small. We set the time step to At = 10~^ GeV~^. 

In two-body final states, one expects no black-hole, 
and hence no black-hole decay by emission of Hawking 
radiation. The generator therefore proceeds directly to 
the final burst phase. 



FINAL BURST 



In the absence of a self-consistent theory of quantum 
gravity, the last stage of the evaporation cannot be de- 
scribed accurately. Once the mass of black-hole becomes 
close to the fundamental scale M», the classical black- 
hole solution can certainly not be used anymore. We 
adopt a scenario in which the final stage of evaporation 
is a burst of particles which conserves energy, momentum 
and all of the gauge quantum numbers. For definiteness, 
we assume the remaining black-hole will decay into the 
lowest number of Standard-Model particles that conserve 
all quantum number, momentum and energy. 

A black-hole with electromagnetic charge Q'''' and 
color- vector c^''' ~ {r^^ ^b^^ ,g^^) will be taken to emit 
Af_i/3 down-type quarks (i.e. d,s or b quarks), 7V2/3 
up-type quarks (u,c, or t), iV_i charged leptons and W 
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bosons , Ngi gluons, and Nn non-charge particles (ic. 7, Z 
an Higgs). We use the following procedure to determines 

Nburst = iN_l/3,N2/3,N.l,Ngi,N^). 

Step 1: preliminary solution: 

• Search all possible solutions with iV„ = 0. 

• Choose the minimum number of particles as pre- 
liminary solution. 

Step 2: Actual charged/colored emitted particle count: 

• The preliminary Nh^rst = 
{N_i/3,N2/3,N_i,Ngi,Nn) having now been 
determined. If the minimum number of solution 
is less than 2, we then add Nn to keep the total 
number equal to 2. Later we choose one of them 
randomly according to the degrees of freedom of 
each particle. 

• After obtaining the number of emitted particles, 
we randomly assign their energies and momenta, 
subject to the constraint that the total energy and 
momentum equal that of the final black-hole state. 
We currently neglect any bulk components of the 
final black-hole momentum. 

VI. INPUT AND OUTPUT 



The input parameters for the generator are read from 
the file parameter.txt, see FigfTil 

1. Number_of ..simulations: sets the total number 
of black-hole events to be simulated; 

2. Center_of _jnass__energy_of __protons: sets the 
center-of-mass energy of the colliding protons in 
GeV; 

3. M__ph: sets the fundamental quantum-gravity scale 
(A/,) in GeV; 

4. Choose_a__case: defines the extra dimension 
model to be simulated: 

(a) 1: non- rotating black-holes on a tensionless 
brane with possibility of fermion splitting, 

(b) 2: non-rotating black- holes on a brane with 
non-zero positive tension, 

(c) 3: rotating black-holes on a tensionless brane 
with d=5, 

(d) 4: two-particle final-state scenario; 

5. number_of __extra_dimensions: sets the number 
of extra dimensions; this must equal 2 for brane 
with tension (Choose_a__case=2); 



6. nuiiiber_of __splitting_dimensions: sets 
the number of extra split-fermion dimensions 
(Choose_a__case=l); 

7. extradimension_size: sets the size of the mini- 
bulk® in units of 1/TeV (Choose_a__case=l); 

8. tension: sets the deficit-angle parameter i? [flfisj 
(Choose_a__case=2); 

9. choose_a__pdf __f ile: defines which of the differ- 
ent CTEQ6 parton-distribution functions (PDF) to 
use; 

10. Minimum_mass: sets the minimum mass M^in in 
GeV of the initial black- holes; 

11. f ix_time_step: If equal to 1, then code uses the 
next parameter to determine the time interval be- 
tween events; if equal to 2 then code tries to opti- 
mize the time step, keeping the probability of emit- 
ting a particle in any given time step below 10%. 

12. time_step: defines the time interval At in GeV~^ 
which the generator will use for the black-hole evo- 
lution; 

13. Mass_loss__f actor: sets the loss factor for the 
energy of the initial black-hole, as defined in equa- 
tion[2l 

14. momentum_loss_f actor: defines the loss factor 
< /p < 1 for the momentum of initial black-holes 
as defined in equation!^ 

15. Angular_momentum_loss_f actor: sets the loss 
factor < /l < 1 for the angular momentum of 
initial black-holes s defined in equation[29l 

16. Seed: sets the seed for the random-number gener- 
ator (9 digit positive integer); 

17. L suppression: chooses the model for suppress- 
ing the accumulation of large black-hole angular 
momenta during the evolution phase of the black- 
holes (cf. discussion surrounding equations [5^115^ ; 

• 1: no suppression; 

• 2: A Area suppression; 

• 3: Jbh suppression; 

• 4: A J suppression; 



This is the distance between fermion branes where only gauge 
bosons and Higgs field can propagate in split-fermion brane sce- 
nario. 
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FIG. 14: Parameter.txt is the input file containing the parameters that one can change. The words in parentheses are the 
parameters that are used in the paper. 
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18. angular_inoineiituin_suppression_f actor: de- 
fines the phenomenological angular-momentum 
suppression factor, (cf. discussion surrounding 
equation [52ll54|) ; 

19. charge_suppression: turns the suppression of ac- 
cumulation of large black-hole electromagnetic and 
color charge during the black-hole evolution pro- 
cess on or off (cf. dicussion surrounding equation 

• 0: charge suppression turned off; 

• 1: charge suppression turned on; 

20. charge_suppression_f actor: sets the electro- 
magnetic charge suppression factor, (q, in 1371 

21. color_suppression_f actor: sets the color 
charge suppression factor, (^3 in 1371 

21-94 (odd entries:) the widths of fermion wave func- 
tions (in AI^^ units); and (even entries:) centers 
of fermion wave functions (in units) in split- 

branc models, represented as 9-dimensional vectors 
(for non-split models, set all entries to 0). 

When the code terminates, the file output.txt with all 
the relevant information (i.e. input parameters, cross sec- 
tion) is output to the working directory. This file contains 
also different segments of information about the gener- 
ation of black-holes which are labelled at the beginning 
of each line with an ID word (Parent, Pbh, trace, Pem, 
Pemc or Blast): 

• Parent: identifies the partons whose collision re- 
sulted in the formation of the initial black-hole (see 

Fig. [m. 

— column 1: identifies the black-hole; 

— column 2: PDGID code of the parton; 

— column 3: energy of the parton; 

— columns 4-6: branc momenta of the parton. 

• Pbh: contains the evolution of the charge, color, 
momentum and energy of the black-holes, and, for 
rotating black- holes, their angular momentum (cf. 
Fig. [TT]). 

— column 1: identifies the black-hole; 

— column 2: time at which the black-hole emit- 
ted a particle; 

— column 3: PDGID code of a black-hole; 

— column 4: three times the electromagnetic 
charge of the black-hole; 

— columns 5 to 7: color-charge vector compo- 
nents of the black-hole; 



— columns 8: energy of the black-hole in the lab- 
oratory frhame; 

— columns 9 to 11: brane components of 
the black-hole momentum in the laboratory 
frame; 

— columns 12 to (8-|-d): bulk components of the 
black-hole momentum; 

— column (9+d): angular momentum of the 
black-hole, in the case of rotating black-holes; 
empty otherwise. 

• trace: contains the evolution history of the black- 
holes' positions (cf. Fig. [T8|) : 

— column 1: identifies the black-hole; 

— column 2: the times at which the black-hole 
emitted a particle; 

— columns 3 to 5are the brane components of the 
black-hole position vector when the black-hole 
emitted a particle; 

— columns 6 to (2-|-d): the bulk components of 
the black-hole position vector, when the black- 
hole emitted a particle. 

• Pem: contains a list of the black-holes, with the 
history of their evolution (cf. Fig. [TO]) : 

— column 1: identifies the black-hole; 

— column 2: the times at which the black-hole 
emitted a particle; 

— column 3: PDGID code of the emitted parti- 
cle; 

— column 4: three times the charge of the emit- 
ted particle; 

— columns 5 to 7: color-vector components of 
the emitted particle; 

— columns 8: energy of the emitted particle in 
the laboratory frame; 

— columns 9 to 11: branc components of the mo- 
mentum of the emitted particle, in the labo- 
ratory frame; 

— columns 12 to (8-|-d): bulk components of the 
momentum of the emitted particle. 

• Pemc: contains the same information as Pem, but 
in the center-of-mass frame of the collision. 

• Elast: contains the same information as Pem for 
the particles emitted in the final decay burst of the 
black-hole. Column 12 and onwards are omitted as 
these particles have no bulk momentum. 
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FIG. 15: Output.txt: There are three parts to this file. The first part is a copy of parameter.txt. The second part includes 
information about the black-hole and the emitted particles. The first column identifies the what type of information each row 
is supplying - parent is information about the two incoming partons; Pbh is information on the energy and momenta of the 
produced black-holes; trace describes the location of the black-holes; Pem characterizes the emitted particles in the lab frame; 
Feme characterizes the emitted particles in the center-of-mass frame; Elast describes the final burst. 
The third part of the file is the black-hole production cross-section as inferred from the events in this generator run. 



VII. RESULTS 

We choose the following parameters for the distribu- 
tions shown in this section and normaUze them to an 
integrated luminosity of lOfb^^, unless otherwise stated. 
The values of the parameters are chosen to be the same 
as in figure I14[ except where a parameter is varied to 
study its effect. 

• Number_of ..simulations = 10000; 

• Center_of _jnass__energy_of __protons = 14000 Ge 
V; 

• M__ph = 1000 GeV; 



• extradimension_size = lOTeV^ ; 

• choose_a__pdf __f ile=0; 

• Minimuni_mass=5000GeV; 

• time_step=10"^ GeV"^ 

• Mass_loss__f actor=0.0; 

• momentuni_loss_f actor=0.0; 

• Angular _momentum_loss_factor=0. 2; 

• Seed=123589341; 

• L suppression=l; 
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FIG. 16: Lines in the output file headed by the ID = Parent contain information about the initial partons which formed the 
black-hole. 
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FIG. 17: Lines in the output file headed by the ID — Pbh contain the energies and momenta of the black-holes for each emission 
step. In case of rotating black-holes, the last column in the line is the angular momentum. 



• charge_suppression=l; 

• charge_suppression_f actor irrelevant since 
charge_suppression=l; 

• color_suppression_f actor irrelevant since 
charge_suppression=l; 

• widths of fermion wave functions = 1 M^^ ; 

• centers of quark wave functions: 
(10-2/3,0,0, ....)(GeV"^) (i.e. all quarks have 
Gaussian wave functions centered on a point 
displaced from the origin by 10^^ /3GeV~^ in the 
first splitting direction); 

• distribution center of leptons: 
(-10-2/3,0, 0,....)(GeV-i). 

In this section, we present some distributions of prop- 
erties of the initial and evolving black-holes and of the 
particles which are emitted by them during the Hawking 
radiation and final-burst phases. 



A. Mass of the Initial Black-Holes 

Figures HU [55] and [52] show the the initial black-hole 
mass distribution for three different extra dimension sce- 
narios: non-rotating black-holes on a tensionless brane, 
non-rotating black-holes on a non-zero tension brane and 
non-rotating black-holes with split fermion branes respec- 
tively. Because we chose 5 TeV as the minimum mass of 
the initial black-hole, the distributions have a cut off at 
5 TeV. 



B. Movement of Black-Holes in the Bulk 

The generator includes recoil of black-holes due to 
Hawking radiation. The recoil modifies the spectrum due 
to the Doppler effect. Even if the effect is small, the high 
energy tail of the emitted particle's energy spectrum is 
longer than for pure Hawking radiation. Fig. [53] shows 
the random motion in the mini-bulk for 10, 000 black- 
holes as consequence of recoil. While most of the black- 
holes remain on the brane where they were formed, a 
significant number of them are capable of drifting all the 
way to the lepton brane. There are also a few events 
where the black holes leave the mini-bulk completely. 
Since Standard-Model charges are confined to the mini- 
bulk, a black-hole needs to carry zero charge in order to 
be able to leave the mini-bulk. Once out of the mini-bulk, 
a black-hole cannot emit Standard-Model particles any- 
more. Models with an additional bulk Z2 symmetry (e.g. 
Randall Sundrum models) do not allow for a black-hole 
recoil from the brane [l^l • Unfortunately, the number of 
black-holes which escape the mini-bulk is so small that 
experimentally we are unlikely to be able to distinguish 
between models on this basis. 



C. Initial Black-Hole Charge Distribution 

Most of the initial black-holes are created by u and d 
quarks. Denote by N^q the number of black- holes that 
have electromagnetic charge Q. Fig. [25] is a histogram 
of 3(5 for Us = 0, = 4, and = 7 in d = 10 space. 
Since at these parton momenta, there are roughly twice 
as many u-quarks in a proton as d-quarks, we expect that 
most of the black- holes have 3Q = 4. i.e. are made of two 
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FIG. 18: Lines in the output file headed by the ID — trace contain the location of the black-hole fo reach emission step. 




FIG. 19: Lines in the output file headed by the ID — Pern contain the types of the emitted particles, their energies and momenta 
in the lab frame and the times of their emission. 




FIG. 20: Lines in the output file headed by the ID — Elast contain the types, energies and momenta of particles of the final 
burst. 
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FIG. 21: Mass distribution of initial black-holes (rotating and 
non-rotating) on a tensionless brane for various numbers of 
extra dimension. 



FIG. 23: Mass distribution of initial (non-rotating) black- 
holes on a tensionless brane for d = 10 and different numbers 
of split-fermion branes. 




FIG. 22: Mass distribution of initial (non-rotating) black- 
holes on a non-zero tension brane for B — 1.0, B — 0.8 and 
B = 0.6. 

u-quarks (fuu)- One docs indeed see a large peak at SQ — 
4 in figure A second peak at 3Q ~ 1 corresponds to 
black-holes made of one d and one u-quark, or from one 
anti-d and one gluon. Since there are only a few gluons 
or anti-quarks at these momenta, j\ ~ jud- 

Similarly, the small peak at 3Q = — 2 is predominantly 
dd and not ti- gluon. 

We exp ect that Jud — 2y/ fuufdd, and thus /i ~ 
2\/ fif-2- This relation is roughly satisfied in Fig. [25l 
In the split-fermion case, since gluons can move in the 
mini-bulk, there is a further suppression of the gluon con- 
tribution due to the wave-function-overlap suppression 
between the gluons and fermions. In particular SQ = 2 
and 2>Q = —1, which are dominated by gluon-quark col- 
lisions, are suppressed, as can be seen in Fig|25l For 
a large number of split dimensions, there are almost no 



FIG. 24: black-hole movement in the mini-bulk due to recoil. 
X\, Yi axe coordinates in two extra dimensions. The red 
circle indicates the width of the quark brane. The blue circle 
indicates the width of the lepton brane. The black lines are 
black-holes traces. The size of the mini-bulk is lOTeV^^ x 
lOTeV^V black-holes with non-zero standard model gauge 
charges bounce back from the wall of the mini-bulk, black- 
holes with zero Standard Model gauge charges can leave the 
mini-bulk. 

gluon-gluon or giuon-quark black- holes. The decline in 
the giuon-quark configurations accounts for the simulta- 
neous rise in the fraction of quark-only configurations (i.e. 
3Q = 4,l,-2). 

D. Initial Black- Hole Color Distribution 

The colliding partons that form the black-hole carry 
gauge charges, in particular color and electromagnetic 
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FIG. 25: Electromagnetic charge distribution of the initial 
black-holes. 



FIG. 26: Initial color distribution of the created black-holes. 
The vertical lines are error bar. 



charge. From the PDFs [S^l we see that, at the relevant 
parton momentum, most of the partons are u and d type 
quarks - essentially the valence quarks. Contributions 
from "the sea" - other quarks, antiquarks, gluons and 
other partons - are subdominant. We therefore estimate 
the distribution of the colors of the initial black-holes to 
be 



7V,(0) : Ar,(l) :iV,(2) = 4:4: 1. 



(55) 



Here Ni (p) is the number of black-holes whose i-th color- 
vector component (i=l is red, i=2 blue, i=3 green) 
has the value p. This agrees very well with the graph 
in Fig. [21 N.i{-1) and N^{-2) refer to black-holes 
created from collisions involving gluons or anti-quarks. 
Their numbers are hard to estimate, but we expect that 
Ni{-2) < /i(-l) < Ni{0 < p < 2), again consistent 
with Fig. UHl 

At energy scales accessible at the LHC, the color dis- 
tribution of black-holes in different brane world models 
does not differ from each other significantly. However, 
were Mmin significantly lower (at or below 1 TeV) , then 
black-hole production by gluon-gluon scattering would 
be more important, significantly altering the color distri- 
bution, and making it more sensitive to fermion brane- 
splitting (which lowers the gluon-gluon contribution). 



E. Evolution of Black-Hole Color and Charge 
during the Hawking Radiation Phase 

Figure [27] shows the color distribution of the black- 
holes which they accumulate during the evaporation 
phase. From equation I55i the expected average initial 
color of the black-holes is 2/3. Since the colors of emit- 
ted particles are assigned randomly, we expect the cumu- 
lative color distribution (CCD) to be symmetric around 
2/3 and peaked at the value. This is indeed what we find. 
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FIG. 27: Cumulative color distribution for non-rotating black- 
holes on a tensionless brane with d = 4 and no fermion brane 
splitting. Histogram with the black squares (open circles) is 
with (without) color suppression. 



The width of the CCD depends on the total number of 
particles emitted by the black-hole during its evaporation 
phase. 

As discussed above, we allow for the possibility of sup- 
pressing particle emission which increases the charge, 
color or angular momentum of the black-hole excessively 
(cf. discussion around equations [571 and . Figure 

\Tl\ shows also the cumulative black-hole color distribu- 
tion where we suppressed the accumulation of large color 
charges during the evaporation phase. In order to am- 
plify the effect of color suppression, we have set /a = 20 
instead of the expected /a ~ 0.1. We see that the num- 
ber of black-holes with a color charge larger than 1 is 
decreased. 
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F. Number of emitted particles 
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FIG. 28: Number of particles emitted by a non-rotating black- 
hole on a tensionless brane prior to the "final burst" as func- 
tion of number of extra dimension. Here Us = 0. The error 
bars denote one standard deviation range. 
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FIG. 29: Number of particles emitted by a non-rotating black- 
hole in the split-fermion model prior to the "final burst" as 
a function of the number of brane-splitting dimensions, with 
d = 10. Error bars denote one standard deviation range. 

Figures [55] through [3T] show the number of particles 
that are emitted by a microscopic black-hole during the 
decay process before its final burst for a variety of models. 

In the single tensionless brane model (Fig. [28]) . the 
number of emitted particles first increases with the num- 
ber of dimensions for a non-rotating black-hole, but then 
decreases. This behaviour is a result of the complicated 
interplay of a number of effects: the horizon size of a 
black-hole of a given mass as a function of d, and its 



i4 



lio 



E 
m 

■5 

!r 

E 



4- 




finite tension brane model 
d=5 n =2 
M =1TeV 

M =5TeV M =14TeV 

vertical line is standard deviation 



04 



0.5 



0.6 



0.7 



— r 

OS 



0.9 



1.0 



b; 



FIG. 30: Number of particles emitted by a non-rotating black- 
hole on a non-zero tension brane prior to the "final burst" as 
function of defizit angle parameter B with d = 5 and Us = 2. 
Error bars denote one standard deviation range. 
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FIG. 31: Number of particles emitted by a rotating black- 
hole prior to the "final burst" as function of number of extra 
dimension with Us = 0. Error bars denote one standard devi- 
ation range. 



effect on the particle emission spectra; the dependence 
of the Hawking temperature on r^/f"^; the existence due 
to energy-momentum conservation of an upper limit of 
]\,fbh^2 on the energy of an emitted particle. The location 
of the peak will shift as a function of the input parameter 
M„iin, the minimum initial mass of a black-hole. 

In the split fcrmion model (Fig. [55)) . the number of 
emitted particles decreases with the number of extra di- 
mensions. This is because, even for a fixed Hawking tem- 
perature, the average energies of emitted gauge bosons 
and scalar fields increase as increases. 

In the model with a finite-tension brane (Fig. 
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FIG. 32: Number of particles emitted at the final burst: the 
average number of final burst is about 3.4. 



the number of particles decreases as the parameter B 
increases, i.e. with decreasing tension. As the tension 
increases, B gets smaller but the horizon radius of the 
black-hole increases. The Hawking temperature therefore 
decreases, and, as a consequence, the average energy of 
emitted particles falls. More particles will therefore be 
emitted in the evolution of the black-hole. 

For rotating black-holes (Fig. [31]) (on a tensionless, 
unsplit brane), the number of emitted particles first in- 
creases, then decreases, and finally reaches a plateau. 
This is due to similar reasons as for non-rotating black- 
holes. Compared to non-rotating black-holes of the same 
mass, rotation shifts the energy of emitted particles to 
higher values because it decreases the horizon radius 
and increases the emission of higher angular-momentum 
modes. This decreases the total number of emitted parti- 
cles. It also means that the effect of the upper kinematic 
limit of M^^ 12 on the emitted particle's energy is mag- 
nified. 

Figure [32] shows the number of particles that are emit- 
ted at the final burst stage for a non-rotating black-hole 
on a tensionless brane with d = 5 and = 0. The 
average number of emitted particles is about 3. During 
the Hawking radiation phase a black-hole emits about 10 
particles, so approximately 30% of the emitted particles 
will be from the final burst stage. In the case exam- 
ined in Figure [321 we did not include suppression of large 
black-hole color or electric charge. Thus some black- holes 
acquire large color and electric charges by the end of the 
Hawking radiation phase. These black-holes then must 
decay into a large number of particles (> 5) in the final 
burst. 



G. Energy Distributions of the Emitted Particles 

Once formed, the black-holes decay by emission of 
Hawking radiation, a process which continues until the 
mass of the black-hole falls to the fundamental quantum- 
gravity scale. At this stage we chose the black-holes to 
burst into a set of Standard-Model particles as described 
in seetion[Vl The observable signatures of the decay will 
depend on the distributions of energy, momentum and 
particle types of the emitted particles. 

Figures [33] through [36] show the relation between the 
mass of the evolving black-hole and the average energy 
of emitted particles for different extra dimension models. 
The error bars denote 1/VN times the standard devia- 
tion of the mean energy. The minimum mass of the initial 
black-hole is taken to be Mmin = 5 TeV. 

We see from figure [551 that, for a single (i.e. unsplit) 
brane, when Mf,/j >> Af*. a black- hole in d = 10 emits 
higher energy particles than a black-hole in lower dimen- 
sions. For black- hole masses closer to M*, the highest 
energy particles are emitted when the dimensionality of 
space is low, i.e. d ~ A. This reversal can be understood 
from figure [H In the LHC energy range, the curves of 
Hawking temperature as a function of black-hole mass 
for different dimensions cross. At high mass, high d ex- 
hibits the highest Hawking temperature; at low mass, 
low d exhibits the highest Hawking temperature. It is 
easy from this figure to estimate the number of emitted 
particles, and to roughly reproduce the results shown in 
figure [281 

The main difference between the curves for different d, 
comes from the changing size of the black- hole horizon. 
For low d, the horizon radius increases more quickly with 
the mass than for higher d, as seen in figure[Tl The Hawk- 
ing temperature of the black-hole is inversely propor- 
tional to r/i . So long as the Hawking temperature remains 
well below the black-hole mass (here for M^h ^ 2TeV), 
the energy of emitted particles decreases as the mass in- 
creases. However as equations (|imi2p show, this change 
is slow for high d. By d = 10, the energy of emitted 
particles is almost constant from 2TeV to hTeV . 

The increase of the energy of emitted particles stops 
at about 2 TeV and then a decrease begins. The rea- 
son, as stated above (equation ([5^ ). is that by energy- 
momentum conservation, a black-hole can only emit par- 
ticles with less than half of its mass. 

In the split-fermion model (Fig. [34]) . the average en- 
ergy of the emitted particles increases as the number of 
dimensions in the mini-bulk increases. The energy shift 
comes from the gauge bosons and scalar fields, which ac- 
cess the higher-dimensional phase space of the mini-bulk 
(cf. appendix A: figure [50] through [61]) . 

For the brane with non-zero tension (Fig. [35]) . the 
radius of the black-hole increases with tension hence the 
energy of emitted particles decreases with tension. 
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FIG. 33: Average energy of the particles emitted by (non- 
rotating) black-holes on unsplit branes versus the mass of the 
black-hole at the time of emission. Here d — 4, d = 7 and 
d = 10. Note that the final burst is not included, because it 
occurs when the mass of the black-hole is less than 1 TeV. 
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FIG. 35: Average energy of the particles emitted by (non- 
rotating) black-holes on a brane with tension, versus the mass 
of the black-hole at the time of emission. Here d = 5, Us = 2 
and B = 1, B = 0.8, B = 0.6. 
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FIG. 34: Average energy of the particles emitted by (non- 
rotating) black-holes on split branes versus the mass of the 
black-hole at the time of emission. Here d = 10 and Us = 0, 
Us — S, Us = 7. 

For a rotating black-hole (Fig. [36]), angular momentum 
decreases the size of the horizon. Thus since black-holes 
are typically formed with some initial angular momen- 
tum, they emit higher energy particles than non-rotating 
black-holes of the same mass. However, the black-hole 
tends to shed its angular momentum rapidly as it emits 
particles. This increases the horizon size, lowers the 
Hawking temperature, and lowers the average energy of 
the emitted particles. The rapid shedding of angular mo- 
mentum thus leads to a drop in the average emitted par- 
ticle energy around A'f,nin- 

If one compares rotating with non-rotating black- holes, 
one finds that the energy of the emitted particles is al- 
ways larger for the rotating black-holes. By the time 



FIG. 36: Average energy of the particles emitted by (possibly 
rotating) black-holes versus the mass of the black-hole at the 
time of emission. Here d = 4, d = 7 and d — 10, with a 
tensionless unsplit brane. 

the mass of the black-hole has dropped well below A/min 
(here to approximately 1 — 2 TeV) , almost all of the angu- 
lar momentum has been lost and the difference between 
the rotating and non-rotating black-holes is small. 

As a Standard-Model particle is emitted from a black- 
hole, this particle may carry some momentum in the ex- 
tra dimensional directions. As an observer in 3-space 
observes this particle, he/she will find the apparent mass 
of the particle is 



Af„fc= ^A/2 + ^F2 (56) 

Here Mob is the observed mass, Mp is the true mass of 
the particle, and Pg is the particles extra-dimensional 
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FIG. 37: Energy distribution of emitted particles in the Hawk- 
ing radiation step for single-brane non-rotating black-hole. 
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FIG. 38: Energy distribution of emitted particles at the final 
burst step for single-brane non-rotating black-hole. 



momentum. Clearly Mob > Mp. A Standard-Model par- 
ticle, however, cannot leave the Standard Model brane. 
Its extra-dimensional momentum must therefore be ab- 
sorbed by the brane or carried away by bulk particles 
(such as gravitons). We therefore calculate an emitted 
particle's "energy on the brane" according to: 
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(57) 



where Pi is the regular 3-momentum. We will assume 
that the shedding of extra-dimensional momentum is 
rapid; and henceforth we will refer to Eb (rather than 
the initial emission energy) as the energy of the emitted 
particle. 

Figure [37] shows the energy distribution of the particles 
from Hawking radiation in the singlc-branc model. The 



cross section for black-hole production increases with d 
(figure[7|) • The area under the curves also increase with d. 
The peaks of the curves are around 200GeV to 400GeV. 
One can compare, for example, the energy distributions 
for 0? = 9 and d = 10. A black-hole in high r d tends to 
emit particles with higher energy, so the curve for c? = 10 
has a longer higher energy tail than for d = 9. 

Figure[551shows the energy distribution of the particles 
from the final burst in the single-brane model. The en- 
ergy these particles share is much smaller than the energy 
in the earlier Hawking radiation phase. The peak in the 
energy distribution of these particle is around 200GeV 
to 300GeV. The tails extend just to ITeV, which is the 
mass at which the black hole is taken to be unstable and 
undergo its final burst. 

Figures [31] through [32] show the energy distribution of 
emitted particles (including final burst particles) in the 
various models that BlackMax can simulate. 

Figure [33] shows the energy distribution in the single- 
brane model. The cross section of black-hole increases 
with d (figure]?]). The area under the curves also increases 
with d. The peaks of the curves are around 200 GeV to 
400 . Again comparing d = 9 and d = 10, a black-hole 
in higher d tends to emit particles with higher energy, so 
the curve for d = 10 has a longer high energy tail than 
for d=9. 

Figure HD] illustrates the split-fermion model. In this 
figure, we keep the total number of dimensions d fixed 
but change the dimensionality Ug of the mini-bulk. This 
affects the spectra of only the gauge boson and scalar 
fields, as only their propagation is affected by the the 
mini-bulk's dimensionality. (Gravitons propagate in the 
full bulk; other Standard-Model particles propagate only 
on the brane.) As explained above, these spectra will 
shift to higher energies as the number of splitting dimen- 
sions is increased. One can see that the curve in Ug = 7 
has the longest high energy tail. 

Figure JJl] illustrates the brane with tension model. 
The energies of the emitted particles shift to lower en- 
ergy as B decreases. 

The energy distribution of emitted particles for rotat- 
ing black-holes (figure has the same general charac- 
teristics as the distribution for non-rotating black-holes. 
Angular momentum causes a black-hole to tend to emit 
higher energy particles than a non-rotating black-hole. 
The curves have longer higher energy tails than the non- 
rotating black-holes. 

Figure 23] shows the Energy distribution of emitted 
particles from two-body final-states scenario. The energy 
of the emitted particles is about the half of the incoming 
partons. 

In figure[44]we show the energy distribution of different 
types of particles in the d = b single-brane model. The 
area of each curve is dependent on the degree of freedom 
of each particle and its power spectrum. One can com- 
pare the ratio of the same type of particles. For example. 
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FIG. 39: Energy density distribution of emitted particles 
from non-rotating black-holes on a tensionless brane with 
no splitted fermion branes. Shown are the distributions for 
4 < d < 10. The spectra include the final burst particles. 



FIG. 41: Energy distribution of emitted particles for non- 
rotating black holes on a non-zero tension brane with B — 1, 
B — 0.8 and B = 0.6. The spectra include the final burst 
particles. 
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FIG. 40: Energy density distribution of emitted particles from 
non-rotating black-holes on a tensionless brane with fermion 
brane splitting. Shown are the distributions for d = 10 and 
< Ha < 7. The spectra include the final burst particles. 



the area of gluons should be 8 times as large as the are 
a of photons. It is roughly the same as what the figure 
shows. 



H. Pseudorapidity Distributions of the Emitted 
Particles 



FiguresBSl to H51 show the pseudorapidity distributions 
of the emitted particles for different extra dimension sce- 
narios. 

Most of black-holes are made of two u quarks and have 
charge 3/4. For the shown figures we did not include 
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FIG. 42: Energy distribution of emitted particles for rotating 
black-holes for 4 < d < 10. 



charge suppression, because of that the black-holes tend 
to emit the same number of particles and antiparticles 
during Hawking radiation phase. This can be seen from 
the curves without final burst in figure H51 through figure 
l48l The majority of the black-holes are positive charged 
and will tend to emit positive particles in the final burst. 
That is why there are more positrons than electrons for 
the distributions which include the final burst particles. 

Figure WS\ shows the pseudorapidity distribution in the 
two-body final-state scenario. The distribution is much 
wider than the equivalent distribution from Hawking ra- 
diation. In this model we do not consider the angular 
momentum of the black hole. We therefore take the de- 
cay process in the two-body final state scenario to be 
isotropic in the ordinary spatial directions, just as in 
other models without angular momentum. In the center 
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FIG. 43: Energy distribution of emitted particles for two- 
body final states. 
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FIG. 45: Pseudorapidity distribution of charge leptons and 
anti leptons for non-rotating black- holes on a tensionless 
brane with and without the final burst particles; d = 5. 



10 -J 




■ Changedlepton 

quark 

■ ■ ■ photon 
-■- WandZ 
- - Hgas 

^uon 

■ neutrino 
d=5 n =0 M=1TeV 

M ._ = 5 TeV M =14TeV 
LLminosity: lOfb ' 



1000 2000 
Enengyof emitted particles (GeVJ 



3000 



FIG. 44: Energy distribution of each particle type in the d = 5 
single brane model. 



of mass frame, particles are therefore emitted in direc- 
tions uncorrelated with the beam direction. Neverthe- 
less, because the threshold energy of the two-body final- 
state model is much lower than that of other models, the 
intermediate state tends to have a higher velocity down 
the beam-pipe. The decay products are therefore emit- 
ted with larger pseudorapidity. In truth, one may expect 
that the intermediate state of the two-body final-state 
scenario has non-negligible angular momentum. How- 
ever, since the intermediate state is not a real black hole, 
it is unclear exactly what role the angular momentum of 
the intermediate state plays. 

In the final state scenario, the momenta of the two 
emitted particles are correlated with the initial parton 
momenta, and hence the pseudorapidity distribution of 
the emitted particles reflects that of the initial partons. 
The ratio of the number of events for pseudorapidity bc- 
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FIG. 46: Pseudorapidity distribution of quarks and anti 
quarks for non-rotating black-holes on a tensionless brane 
with and without the final burst particles; d — 5. 



tween and 0.5 divided by that for pseudorapidity be- 
tween 0.5 and 1 is about 1.1. This is much higher than 
the asymptotic QCD value of 0.6, as predicted by [3, [2^. 
If the ratio is found not to equal 0.6, then this would sug- 
gest new physics beyond the Standard Model. 



I. Emitted Particle Types 

Table Hill shows, for a variety of representative extra- 
dimension scenarios the fraction of emitted particles 
which are of each possible type - quarks, gluons, 
(charged) leptons, (weak) gauge bosons, neutrinos, gravi- 
tons, Higgs bosons and photons. One notable feature is 
that the intensity of gravitons relative to other particles 
increases with the number of extra dimensions. Note that 
the absence of gravitons in the case of a rotating black- 
hole is not physical, but rather reflects our ignorance of 
the correct gray-body factor. 
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FIG. 47: Fseudorapidity distribution of charged leptons and 
anti-leptons for rotating black-holes, on a tensionless brane, 
with and without the final burst particles. Here d = 5. 



FIG. 49: Fseudorapidity distribution for the two-body final- 
state scenario. The distribution in 2-body final states is much 
flatter near 77 = than emission due to Hawking radiation. 
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FIG. 48: Fseudorapidity distribution of quarks and anti- 
quarks for rotating black-holes, on a tensionless brane, with 
and without the final burst particles. Here d = 5. 

VIII. CONCLUSION 

Hitherto, black-hole generators for the large-extra- 
dimension searches at the LHC have made many simpli- 
fying assumptions regarding the model of both our three- 
dimensional space and the extra-dimensional space, and 
simplifying assumptions regarding the properties of the 
black-holes that are produced. In this paper we have 
discussed a new generator for black- holes at the LHC, 
BlackMax, which removes many of these assumptions. 
With regard to the extra-dimensional model it allows for 
brane tension, and brane splitting. With regard to the 
black-hole, it allows for black-hole rotation, charge (both 
electro-magnetic and color) and bulk recoil. It also in- 
troduces the possibility of a two-body final state that is 
not a black-hole. 

Although BlackMax represents a major step forward, 
there remain important deficiencies that will need to be 



addressed in the future. BlackMax continues to insist on 
a flat geometry for the bulk space, whereas there is con- 
siderable interest in a warped geometry Q or in a com- 
pact hyperbolic geometry While BlackMax allows 
for black-hole rotation, the absence of either an analytic 
or a numerical gray-body factor for the graviton in more 
than three space dimensions for rotating black-holes is 
a serious shortcoming that can be expected to materi- 
ally change the signature of black-hole decay for rotating 
black-holes. Other issues include how to properly ac- 
count for the likely suppression of decays that cause a 
black-hole to acquire very "large" color, charge or angu- 
lar momentum. (As oppposed to the somewhat contrived 
phenomenological approach currently taken.) These are 
but a few of the fundamental issues that remain to be 
clarified. 

Despite these (and no doubt other) shortcomings, we 
expect that BlackMax will allow for a much improved un- 
derstanding of the signatures of black-holes at the LHC. 
Work in progress focuses on using BlackMax to explore 
the consequences for the ATLAS experiment of more re- 
alistic black-hole and extra-dimension models. 
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Appendix A: Particle Emission Spectra 

Nowadays gray-body factors can be found in many pa- 
pers. We collect the relevant papers in tabic [D We 
follow these papers to calculate energy power spectra 
in our database except for the split-fermions model (al- 
though we independently confirm the results as well) . We 
perform an original calculation of the spectra of gauge 
bosons and scalar fields in the split-fcrmions model as 
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a function of the number of dimensions in which the 
fermion branes are split. These are shown in figures [50l 
toiB 

Anyone who has a better spectrum can upgrade our 
database. For example, we only calculate the spectra up 
to the I = 9 mode. 
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FIG. 52: Spectra of scalar fields in d=7 space. 



FIG. 50: Spectra of scalar fields in d=5 space. 
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FIG. 53: Spectra of scalar fields in d=8 space. 
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FIG. 51: Spectra of scalar fields in d=6 space. 




FIG. 54: Spectra of scalar fields in d=9 space. 
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FIG. 55: Spectra of scalar fields in d=10 space. 



FIG. 58; Spectra of gauge bosons in d=7 space. 
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FIG. 56; Spectra of gauge bosons in d=5 space. 



FIG. 59; Spectra of gauge bosons in d=8 space. 
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FIG. 57; Spectra of gauge bosons in d=6 space. 



FIG. 60; Spectra of gauge bosons in d=9 space. 




FIG. 61: Spectra of gauge bosons in d=10 space. 
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TABLE I: Literature sources for particle emission spectra 





Tvnp nf blacW Vinlp 


RrflTip TTindpl 




Standard-Model particles 


non-rotating 


unsplit; tensionless 


[18J[19| 


gravitons 


non-rotating 


split / unsplit ; tensionless 


[17| 


Standard-Model particles 


non-rotating 


split /unsplit; with tension 




gravitons 


non-rotating 


split/unsplit; with tension 


[15] 


scalars and gauge bosons 


non-rotating 


split; tensionless 


figures [501I6T] 


fermions 


rotating 


unsplit; tensionless 


[23] [24] 


gauge bosons 


rotating 


unsplit; tensionless 




scalar fields 


rotating 


unsplit; tensionless 


[18] [20j[2]J[2^ 



TABLE IL Degrees of freedom of Standard-Model particles which are emitted from a black hole. For gravitons, the table 
shows 1, because the appropriate growth in the number of degrees of freedom is included explicitly in the graviton emission 
spectrum. Us is the number of extra dimensions in which vector and scalar fields can propagate. 



particle type 


do 


di/2 


di 


d2 


Quarks 





6 








Charged leptons 





2 








Neutrinos 





2 








Photons or gluons 








2 + Us 





Z" 


1 





2 + Us 





W+ and W~ 


2 





2{2 + n,) 





Higgs boson 


1 











Graviton 











1 
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TABLE III: The fraction of emitted particles of different types (including final burst particles) in a variety of extra dimension 
scenarios. *Note that the absence of gravitons in the case of a rotating black hole is due exclusively to our current ignorance 
of the correct gray-body factor. 



Scenario 


quarks 


gluons 


leptons 


gauge bosons 


neutrinos 


gravitons 


Higgs bosons 


photons 


d = 4 lis = non-rotating black hole 


68.21 


10.79 


9.45 


5.72 


3.87 


2.00e-01 


8.99e-01 


8.61e-01 


d = 5 Us = non-rotating black hole 


65.37 


13.29 


9.04 


6.12 


3.76 


4.60e-01 


8.26e-01 


1.13 


d = 6 Us = non-rotating black hole 


63.63 


14.51 


8.93 


6.58 


3.51 


7.59e-01 


7.76e-01 


1.30 


d = 7 Us = non-rotating black hole 


61.25 


15.94 


8.75 


7.17 


3.45 


1.20 


7.89e-01 


1.44 


d = 8 Us = non-rotating black hole 


60.99 


15.94 


8.56 


6.99 


3.35 


1.93 


7.63e-01 


1.47 


d = 9 Us = non-rotating black hole 


59.40 


16.26 


8.26 


7.08 


3.26 


3.48 


7.45e-01 


1.51 


d = 10 = non-rotating black hole 


57.56 


16.15 


7.68 


6.82 


3.17 


6.46 


6.97e-01 


1.46 


d = 10 = 1 split fermions model 


61.58 


19.76 


1.64 


7.30 


3.83e-01 


6.95 


5.88e-01 


1.80 


d = 10 ns = 2 split fermions model 


61.28 


20.40 


1.66 


7.09 


3.49O-01 


6.91 


4.56e-01 


1.86 


d = 10 = 3 split fermions model 


62.10 


20.33 


1.65 


6.76 


4.46e-01 


6.43 


4.15e-01 


1.87 


d = 10 ns = 4 split fermions model 


62.70 


19.76 


1.72 


6.60 


4.66e-01 


6.62 


3.28e-01 


1.80 


d = 10 = 5 split fermions model 


63.67 


19.34 


1.73 


6.14 


5.17O-01 


6.50 


2.65e-01 


1.83 


d = 10 ns = 6 split fermions model 


64.96 


18.24 


1.82 


5.85 


5.56O-01 


6.69 


2.35e-01 


1.65 


d = 10 ns = 7 split fermions model 


66.38 


17.23 


1.83 


5.44 


5.50O-01 


6.68 


2.51e-01 


1.64 


d = 2 Us = 2 B — 1.0 tension brane model 


84.13 


9.39 


1.66 


2.80 


2.02C-01 


9.23e-01 


2.51e-01 


6.30e-01 


d = 2 Us = 2 B = 0.9 tension brane model 


85.33 


8.88 


1.65 


2.49 


1.63e-01 


7.55e-01 


1.90e-01 


5.47e-01 


d = 2 Us = 2 B = 0.8 tension brane model 


86.42 


8.32 


1.54 


2.32 


1.79e-01 


5.70e-01 


1.99e-01 


4.49e-01 


d = 2 Us = 2 B — 0.7 tension brane model 


87.06 


7.93 


1.63 


2.08 


1.75O-01 


5.19e-01 


1.98e-01 


4.00e-01 


d = 2 Us = 2 B — 0.6 tension brane model 


87.92 


7.39 


1.60 


1.94 


1.59C-01 


4.46e-01 


1.88e-01 


3.59C-01 


d = 2 Ua = 2 B = 0.5 tension brane model 


88.51 


7.11 


1.46 


1.84 


1.33e-01 


3.96e-01 


1.86e-01 


3.67e-01 


d = 2 Us = 2 B = 0.4: tension brane model 


89.16 


6.58 


1.58 


1.72 


1.49e-01 


3.54e-01 


1.66e-01 


2.90e-01 


d = 4 ris = rotating black hole 


64.82 


15.41 


7.94 


6.25 


3.50 


0.00* 


6.06e-01 


1.48 


d = 5 ns = rotating black hole 


61.38 


18.11 


7.89 


7.06 


3.30 


0.00* 


5.56e-01 


1.70 


d = 6 Us = rotating black hole 


59.21 


20.38 


7.27 


7.59 


3.08 


0.00* 


5.85e-01 


1.89 


d = 7 Us = rotating black hole 


57.52 


21.82 


7.08 


8.08 


2.87 


0.00* 


5.77e-01 


2.05 


d = 8 71s = rotating black hole 


54.41 


24.02 


6.75 


9.09 


2.77 


0.00* 


6.07e-01 


2.36 


d = 9 Us = rotating black hole 


51.88 


24.05 


7.98 


9.27 


3.55 


0.00* 


7.15e-01 


2.55 


d = 10 ns = rotating black hole 


52.43 


25.67 


6.45 


9.62 


2.52 


0.00* 


6.79e-01 


2.63 


d = 4 ris = two-body final states 


58.88 


16.89 


1.38 


6.08 


11.14 


O.OOe-fOO 


1.91 


3.71 


d = 5 ?is = two-body final states 


57.71 


17.46 


1.37 


6.22 


11.53 


O.OOe-KOO 


1.91 


3.80 


d = 6 ris = two-body final states 


57.16 


17.72 


1.39 


6.27 


11.59 


O.OOe-fOO 


1.93 


3.94 


d = 7 Us = two-body final states 


57.02 


17.91 


1.37 


6.22 


11.59 


O.OOe+00 


1.98 


3.91 


d = 8 ris = two-body final states 


56.58 


18.08 


1.35 


6.31 


11.79 


O.OOe-fOO 


1.92 


3.97 


d = 9 71s = two-body final states 


56.51 


18.02 


1.41 


6.40 


11.75 


O.OOo+OO 


1.96 


3.95 


d = 10 ris = two-body final states 


56.66 


18.06 


1.43 


6.28 


11.74 


O.OOe+00 


1.92 


3.90 



